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Many models of electroweak symmetry breaking have an additional light pseudoscalar. If the 
Higgs boson can decay to a new pseudoscalar, LEP searches for the Higgs can be significantly 
altered and the Higgs can be as light as 86 GeV. Discovering the Higgs boson in these models is 
challenging when the pseudoscalar is lighter than 10 GeV because it decays dominantly into tau 
leptons. In this paper, we discuss discovering the Higgs in a subdominant decay mode where one 
of the pseudoscalars decays to a pair of muons. This search allows for potential discovery of a 
cascade-decaying Higgs boson with the complete Tevatron data set or early data at the LHC. 



I. INTRODUCTION 

The last unexplored frontier of the Standard Model 
is electroweak symmetry breaking, the process by which 
the Higgs field obtains a vacuum expectation value and 
gives mass to the and gauge bosons. One of the 
major goals of current colliders is to discover the Higgs 
boson and understand the dynamics that give rise to 
electroweak symmetry breaking. There have been direct 
and indirect searches for the Standard Model (SM) Higgs 
at LEP and the Tevatron. The current lower bound on 
the Higgs mass, 

ruho > IMA GeV (95% confidence), 

comes from searches at LEP for e^e~ Z^h^ , with the 
SM Higgs decaying to a pair of taus or bottom quarks 
[1]. Recently, combined Higgs searches from the CDF 
and experiments at the Tevatron excluded a SM- 
like Higgs of 169 GeV < m^o < 171 GeV 0. 

While direct searches for the Higgs point towards 
a heavy mass, indirect bounds from electroweak con- 
straints place a limit on how heavy the mass can be. 
In particular, the best fit for a SM Higgs mass is 77 
GeV with a 95% upper bound of 167 GeV [3|. This 
limit comes from measurements of electroweak param- 
eters that depend logarithmically on the Higgs mass 
through radiative corrections. There is tension between 
the direct and indirect measurements; only a narrow win- 
dow of masses for the SM Higgs satisfies both results. 

On the theoretical side, a light Higgs is preferred 
within the Minimal Supersymmetric Standard Model 
(MSSM). Requiring a natural theory and minimizing 
fine tuning drives the Higgs mass below the LEP direct 
bound. In the MSSM, there are two new Higgs chiral su- 
perfields, Hu and Hd, that result in two CP-even scalars 
and , the CP-odd scalar , and the charged Higgs 

after electroweak symmetry breaking. Typically, the 



/i" has Standard Model-like couplings. At the one-loop 
level, the Higgs boson mass is 

m\o — m\a cos^ 2f} 

where at is the dimensionless trilinear coupling between 
the Higgs and top squarks 

At-^lC0t(3 
at = , (1) 

For a moderate 04 ^ 1 and top squarks lighter than 1 
TeV, the Higgs mass is less than 120 GeV [HE]. By 
taking at to "maximal mixing," where the contribution 
from the yl-terms gives the largest contribution to the 
Higgs mass, the Higgs can be as heavy as 130 GeV while 
keeping the top squarks under 1 TeV. Two-loop correc- 
tions can raise the Higgs mass by an additional < 6 GeV 

To avoid fine tuning, the top squarks should not be 
significantly heavier than the Higgs. Even with masses 
at 1 TeV, the Higgs potential is tuned at the few percent 
level. If the top squarks are at 400 GeV, the fine tuning 
of the Higgs potential drops substantially; however, the 
upper limit on the Higgs mass falls to 120 GeV even with 
maximal top squark mixing [B] . This has motivated stud- 
ies giving the Higgs quartic coupling additional contri- 
butions inside the supersymmetric Standard Model [7], 
which usually leads to a less minimal Higgs sector such as 
in the next-to-minimal supersymmetric Standard Model 
(NMSSM). 

Alternate models of electroweak symmetry break- 
ing that can have naturally light Higgs bosons are moti- 
vated by the indirect bounds coming from electroweak 
constraints and the desire to minimize fine tuning in 
the Higgs sector. These models, which often have more 



elaborate Higgs potentials with additional scalar fields, 
allow light Higgs masses, while simultaneously evading 
the LEP direct bound. Frequently, these less mini- 
mal models of electroweak symmetry breaking have ap- 
proximate global symmetries and light pseudo-Goldstone 
bosons that can alter the phenomenology of the Higgs. 
These light pseudo-Goldstone bosons can evade all exist- 
ing limits because they couple very weakly to light flavor 
fermions. 

In this paper, we will focus on such non- minimal 
models of electroweak supersymmetry breaking. We be- 
gin in Section II with a brief discussion of Higgs mod- 
els that contain light pseudo-Goldstone bosons, focusing 
primarily on current experimental constraints. In Sec- 
tion HI, we propose a new search for Higgs bosons that 
cascade decay to pseudoscalars with masses below 10 
GeV. In particular, we find that the Higgs can be dis- 
covered in a subdominant decay mode where one pseu- 
doscalar decays to muons and the other to taus. We 
conclude with a discussion of the expected sensitivity 
to the Higgs production cross section at the Tevatron 
and LHC. The proposed search would allow possible dis- 
covery of a cascade-decaying Higgs with the complete 
Tevatron data set or early data at the LHC. 



II. LIGHT a° MODIFICATIONS TO HIGGS 
PHENOMENOLOGY 

In this section, we explore the couplings of a light 
pseudo-Goldstone boson, or "axion," to the Standard 
Model. We will describe how to analyze a general the- 
ory and we find constraints on the maximal width for 
a Higgs decaying into a light axion. Finally, we will 
analyze how CLEO limits on the direct coupling be- 
tween pseudo-Goldstones and Standard Model fermions 
set constraints on the Higgs width into axions. For Type 
II two Higgs doublet models (i.e., MSSM), the CLEO 
and LEP results place important limits on the light ax- 
ion scenario. 

The tension between the LEP limit on the Higgs 
mass and fine tuning can be reduced if the Higgs branch- 
ing fractions are altered from those of the Standard 
Model [g; . The LEP bound of 114 GeV only applies when 
the Higgs decays dominantly to a bb or t~^t~ pair. While 
adding new, invisible decay modes does not help because 
bounds for such processes are just as strong [5], other 
nonstandard decays remain open possibilities. Consider 
the case where the Higgs decays dominantly to two new 
scalars 0, which in turn decay to SM particles: 

h° ^ (j}<j)-^ {XX){XX). (2) 

For an with SM-like production cross section, this 
process is excluded for m/jo < 110 GeV and X = b [101 
[TTJ[T2]. However, there is an 82 GeV model-independent 



bound from LEP [13j and when X = g,c, t, there are no 
limits for Higgs masses above 86 GeV [111. 

The decay width of a 100 GeV Higgs into Standard 
Model particles is T{h^ SM) ~ 2.6 MeV; because 
the decay width is so small, the Standard Model de- 
cay mode is easily suppressed by the presence of new 
decay modes. Any new light particle with 0(1) cou- 
pling to the Higgs will swamp the decay modes into 
SM particles. Many theories, such as little Higgs mod- 
els and non-supersymmetric two Higgs doublet models, 
have light neutral states for the Higgs to decay into. This 
phenomenon arises when there is an approximate sym- 
metry of the Higgs potential that is explicitly broken 
by a small term in the potential. There is a resulting 
light pseudo-Goldstone boson that couples significantly 
to the Higgs boson. The Peccei-Quinn symmetry of a 
two Higgs doublet model is one such example. In the 
MSSM, it is possible to have a light A°, even with ra- 
diative corrections included More often, there is 
an additional singlet, 5, and an approximate symmetry 
that acts upon the Higgs boson doublets as Hi e*^'' Hi 
with the singlet compensating by S* — > e^^'^^S. During 
electroweak symmetry breaking, S also acquires a vev, 
spontaneously breaking the symmetry; the phase of S 
becomes a pseudo-Goldstone boson and has small inter- 
actions with the Standard Model when {S) » v. 

For specificity, let us consider a two Higgs dou- 
blet model with an additional complex singlet.^ All 
three scalar fields acquire vacuum expectation values: 
Vu — vsinP, Vd — vcos/3, and (S). The interactions 
of the pseudo-Goldstone bosons can be described in the 
exponential basis 

/ cc)+ sin /3 + /i^ cos /3 \ j^i^ 

\ —u) cos p + n sinp J 
S = -L^{S) + s')e'fi>. (3) 

The pseudoscalar fields , ad and get interac- 
tions either through derivative couplings from the kinetic 
terms or through explicit symmetry breaking. One linear 
combination of the pseudoscalars, 

LOza — — a„ sin f3 + Ud cos /?, 

becomes the longitudinal component of the Z*^. The 
two other combinations are physical fluctuations that get 



^ A similar analysis was performed for a one Higgs doublet and a 
complex singlet in |16| . A one Higgs doublet model with a light 
pseudoscalar typically does not have a large branching ratio of 
the Higgs into pseudoscalars. 



2 



mass through symmetry breaking effects in the Higgs po- 
tential. Any terms in the Higgs potential proportional 
to or \S\'^ do not affect the mass or inter- 

actions of the pseudo-Goldstones and there are only a 
handful of possibilities for explicit symmetry breaking. 
As an example, consider adding to the potential a size- 
able coupling 



Vi = XiS^hIhI + h.c. 



(4) 



This will give a weak-scale mass to the following linear 
combination of pseudo-Goldstones: 

— cos Oa{au cos (3 + ad sin (3) — sin 9a- (5) 

The singlet mixing component is given by 

V 

tan 9a = sin 2/3. (6) 

The remaining linearly independent pseudo-Goldstone 
will be massless until the final symmetry is broken. This 
linear combination is 

a° = sin 6'a(a„ cos P + ad sin /3) + cos 9a (7) 

and gets a mass through potentials such as 

V2 = \2S^H^Hd + h.c. (8) 

There will be mixing between a" and A'^ . However, when 
Ai ^ A2, A° and a'^ are nearly mass eigenstates with 
residual mixing proportional to m^o / rri^o . 

It is worth noting that as (5*) <C v, the light pseu- 
doscalar becomes the Peccei-Quinn pseudoscalar and its 
couplings are independent of (S). This particular ex- 
ample has the same symmetry structure as the NMSSM 
near the R-symmetric limit when SH^Hd A-term dom- 
inates the A-term. We will couch our discussions in 
terms of the R-symmetric NMSSM for comparison with 
the literature, but other realizations of the symmetry 
breaking are just as applicable. 

To evade limits from LEP, the Higgs needs a sig- 
nificant branching rate into the pseudoscalar (Fig. 1), 
and one might worry that radiative corrections from the 
Higgs-pseudoscalar interaction might induce a large ra- 
diative correction to the pseudoscalar mass. However, 
the interaction that leads to the Higgs decay into pseu- 
doscalars can occur even if the axion is an exact Gold- 
stone boson through the coupling 



9 

dh 

V 



(9) 



The first interaction preserves the a° 



- e shift sym- 



metry, where Ch is an 0(1) constant and v = 246 GeV 
is the electroweak scale. Because this coupling exists in 
the symmetry-preserving limit, Ch can only depend on 
the vevs of the Higgs fields, the particular charges of the 
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FIG. 1; The branching fraction of the Higgs into pseu- 
doscalars as a function of {S) / sin 2/3 for m/jO = 100 GeV 
when dh = and 1 (solid and dashed lines, respectively). 
The inset shows the minimum value of the branching rate 
into the Standard Model as a function of 771^0 . 



approximate U{1) symmetry, and on the alignment of 
the physical Higgs boson relative to the Higgs vev direc- 
tion. When the physical Higgs boson is in the direction 
of the Higgs vev, /i° = hu sin (3 + hd cos /3, the example 
above gives 



Ch 



: sin^ 9, 



sin^ 2/3 

I v"^ sin^ 2/3 



tan^/3 



(10) 



The second interaction breaks the shift symmetry and is 
proportional to m^o • This term depends on the symme- 
try breaking that gives the axion a mass and is therefore 
model-dependent. When the physical Higgs boson aligns 
with the Higgs vev, the symmetry breaking coupling sim- 
plifies to 



dh = 1 



(11) 



for the potential in Eq. |8] For a symmetry breaking 
potential 



h.c. 



(12) 



dh would be small, arising from the residual mixing be- 
tween and . 

It is possible to increase dh by having multiple terms 
in the potential contribute to m^o , with the pseudoscalar 
mass being less than either of the contributions. For 
instance, with V2 and V2' 



1 



2A2/ sin 2/3 
A2 tan2 e„ 



> 1 if A, 



A2 tan^ 9, 



2 sin 2/3 



-. (13) 



Of course, this is the technical definition of fine tuning 
and when d^, ^ 1, a° has been fine tuned to be light. 
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FIG. 2: Values of (S) / sin 2(3 (GeV) that have been excluded 
through LEP2's search for a Standard Model Higgs [T]. The 
region below m^o = 86 GeV is entirely excluded by the /i" 
4r search [14] , 



We will discount this possibility in our discussion on the 
expected sizes of couplings in this class of theories. 

The partial width of the Higgs into pseudoscalars 
from these interactions is 




quarks. For Type II two Higgs doublet models,^ the 
axion coupling to fermions is given by 



Ant = Wf Jlbja , 



(17) 



where 



9f = sin 9a 




(up-type quarks) 
(down-type quarks/leptons) 



(18) 



(see also [THllSn])- The coupling to up- type quarks is sup- 
pressed by two powers of tan/?. This means that, above 
the 6-quark threshold, the axion will preferentially decay 
to 5-quarks. Below this threshold, it will preferentially 
decay to tau leptons, rather than charm quarks. 

The CLEO bound on the branching fraction of T 
sets a bound on g^, which can be used to set limits on 
the allowed range of gd vs. rrigO (Fig. [s]). The CLEO 
results place the strongest constraints on small values of 
(S) I sin 2/3. The value of the singlet vev is a measure of 
the fine tuning of the theory because it induces a mass 
for the Higgs bosons, m^^ = ^{S)"^ in the scalar potential 
and with an 0(1) coupling, the singlet vev should be less 
than a few TeV to avoid large fine tuning [21 . There 
is some tension, then, between keeping the coupling to 
fermions small and keeping the coupling to the Higgs 



The symmetry-preserving interaction dominates when 



Ch 

2dh. 



m„o 



< 1.5 TeV. 



(15) 



When (S) is less than 1 TeV, the Higgs boson has an 
appreciable width into pseudoscalars (Fig. 1). This is 
precisely the region we are interested in, so we will set 
dfi = for the rest of this discussion. 

Figure [2] shows the values of (S) / sin 2/3 that are 
necessary to evade LEP2's search for a Standard Model 
Higgs [1]. For Higgs boson masses less than 100 GeV, 
the a" has to be fairly strongly coupled to the Higgs 
boson, requiring (S) to be small, which increases the size 
of the coupling to Standard Model fermions. However, 
bounds from the recent CLEO results [17] are strongest 
for large fermion couplings. The CLEO bounds place a 
90% C.L. upper limit on Br(T ja°) Br(a° t+t~). 
For TTigO between 3.5 GeV and 9 GeV, this limit ranges 
from ^ 10~^ — 10""* for the tau decays. The branching 
fraction for radiative T decays is [TH] 



Br(T ^ a"7) 
Br(T-^/.+Ai-) 



F 



(16) 



where F is a QCD correction factor ~ 0.5, Br(T(ls) 
/i^/z^) = 2.5%, and gd is the axion coupling to down 
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FIG. 3: Region of m^o —gd parameter space that has been ex- 
cluded by CLEO to 90% C.L 17. The dashed hues indicate 
values of (S) / sin 2/3 for tan/3 = 2. The shaded region shows 
the minimum values of gd allowed by LEP for an 87-110 GeV 
Higgs. 



^ This case applies to the MSSM and its extensions. For Type I 
two Higgs doublet models, where all Standard Model fermions 
only couple to one Higgs doublet, there is no asymmetry between 
up and down-type quarks in the coupling to axions and this 
typically results in a cot P suppression in the coupling to axions. 
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boson sufficiently large to evade LEP limits without fine 
tuning the a° to be light. 

ft is also possible for LEP to have directly produced 
the Higgs through e+e" Z° ^ h^aP 12J. The LEP 
searches place bounds on the product of the squared 
Z^hPa^ coupling and the branching ratio of the Higgs 
into a Standard Model fermion /: 



sin^ Oa sin^ 2/3 



1 



12 " 



■Br(/i° ^ //)sM < 



nib 



■ sm 



rriho 



(19) 



where Br(/i° //)sm is the Standard Model's branch- 
ing ratio to fermion pairs. There were searches for the 
(6B)(T+r~) final state at LEP, but there were no limits 
for 75 GeV < m^o < 125 GeV. For 125 GeV < m^o < 
165 GeV, the limits were < 0(0.4) which is automat- 
ically satisfied in these models. There were additional 
searches for the {t'^t^){t^t^) final state, but these con- 
straints are even weaker because the Higgs branching 
fraction into taus is a factor of ten smaller than that 
into bottoms. LEP did search for e+e" ^ /i^a^ — > 

(a°a°)a'', but the search for the 6r final state was only 
performed at LEPl and was thus not sensitive to Higgs 
masses above 75 GeV. 



III. h' 







AT HADRON COLLIDERS 



We will now discuss how the Higgs can be dis- 
covered if it decays into a light pseudoscalar a", when 
2mT- < niao ^ 2mf,. In this range, the a" decays pre- 
dominantly into taus and the signature of the Higgs is 
the appearance of 4r events. All existing searches for this 
decay channel have focused on the scenario where two or 
more taus decay leptonically [551 1131 HI] • Currently, the 
ATLAS collaboration is exploring the AfiSv channel and 
CMS is analyzing if^'^T^){fJ-'^T^) [IS]- There are spe- 
cific challenges to the 4t decay channel, however. The 
branching fraction of the taus to leptons is only 33% and 
the pt spectrum of the events is soft because the visible 
lepton carries less than half the momentum of the tau. 
Additionally, it is challenging to reconstruct the Higgs 
and pseudoscalar masses from the final decay products. 



A. Signal 

The primary innovation of the search proposed in 
this paper is to use the subdominant decay of the a'^ into 
two muons, which exists because the a° couples to the 
Standard Model by mixing through the CP-odd Higgs. 
The relative branching ratio for the a" into muons versus 
taus is 



r(aO ^ T+T-) 



(20) 





Signal Efficiency 


Selection Criteria 


Relative 


Cumulative 


Pre-Selection Criteria 


26% 


26% 


Jet veto 


99% 


26% 


Muon ISO & tracking 


~ 50% 


13% 


M^^ < 10 GeV 


98% 


13% 


pi^" > 40 GeV 


76% 


9.8% 


> 30 GeV 


29% 


2.8% 


A<;/>(/i, fir) > 140° 


73% 


2.1% 


AR(p,^i) >0.26 


63% 


1.8% 



TABLE I: Relative and cumulative signal efficiencies due 
to the specified selection criteria. The signal point is a 100 
GeV Higgs decaying to a 7 GeV a" at the LHC. The pre- 
selection criteria include finding a pair of oppositely-signed 
muons, each with I?;] < 2 and pr > 10 GeV. 



The cross section of ft," — ^ 2/i2T depends upon the fol- 
lowing product of branching ratios: 



Ef^r = 2 Br(a° 



/i+/i") Br(a° 



-)• 



(21) 



For tan /3 ^ 4, a 7 GeV pseudoscalar has a 0.4% branch- 
ing ratio into muons and 98% ratio into taus. As m^jO 
goes from the bottom threshold to the tau threshold, e^,- 
varies from 0.8% to 1.5%. The remaining events go into 
hadrons and are divided between the charm and glue- 
glue decay channels. For tan/3 = 2 the branching ratio 
to charms becomes 15% and the branching ratio to taus 
and muons is reduced to 83% and 0.3%, respectively, 
causing e^^ to fall to 0.2% for a 7 GeV a° and 0.5% 
for an a" just above the tau threshold. The events that 
go into hadrons do not typically have significant missing 
energy and do not pass the missing energy cuts. 

Due to the pseudoscalar's small branching fraction 
into muons, this decay channel has not been explored. 
However, the small branching fraction into muons need 
not be a deterrent. The main contribution to the cross 
section for light (~ 100 GeV) SM-like Higgses comes 
from gluon-gluon fusion and can be as high as 2 pb at the 
Tevatron or 50 pb at the LHC [4 . As a result, it is still 
possible to get 300 events with 20 fb~^ at the Tevatron 
(combined and CDF) and 250 events per experiment 
at the LHC (at ^/s = 14 TeV) with 500 ph'^ luminos- 
ity, despite the small branching fraction to muons. We 
ultimately find 0(2%) cumulative efficiency for the sig- 
nal (Table I) , resulting in 95% exclusion limits in certain 
mass windows at the Tevatron. At the LHC, there is the 
possibility for discovery within the first year of running. 

When the Higgs boson decays to two light CP-odd 
scalars a°, the pseudoscalars are highly boosted and 
back-to-back in the center-of-mass frame (Fig. |4]). We 
consider the case where there is a nearly-coUinear pair 
of oppositely-signed muons on one side of the event and 
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FIG. 4: Schematic of Higgs decay chain. The muons and taus 
will be highly boosted and nearly coUinear. It is likely that 
the taus will be reconstructed as one jet. Most of the fir in 
the event will be in the direction of this jet. 



a nearly-coUinear pair of taus on the other, which we 
refer to as a ditau (dir). Each tau has a 66% hadronic 
branching fraction; consequently, there is a 44% proba- 
bihty that both taus will decay into pions and neutrinos, 
which the detector will see as jets and missing energy. 
Even if the taus do not both decay hadronically, there is 
still missing energy, as well as a jet and a lepton, except 
when both taus decay to muons, which occurs ~ 3% of 
the time. The signal of interest is 

pp fJ.~^^J.~ + dir + ]^T, 



where the missing energy comes from the boosted neu- 
trinos and points in the direction of the ditau. Because 
the taus are nearly coUinear, the ditaus are often not 
resolved, leading to a single jet-like object. 

Signal events for a 7 GeV pseudoscalar decaying 
into 2^2r (e^,- = 0.8%) were generated, showered, and 
hadronized using PYTHIA 6.4 Unlike at LEP, the 

overall magnitude of the Standard Model Higgs produc- 
tion cross section is sensitive to physics beyond the Stan- 
dard Model and it is possible to increase the cross section 
by an order of magnitude by adding new colored parti- 
cles that couple to electroweak symmetry breaking. In 
this study, the NNLO Standard Model production cross 
section was used as the benchmark value |1] . 

PGS [27] was used as the detector simulator. Be- 
cause the muons are adjacent, standard isolation cannot 
be used. The muon isolation criteria must be modified 
to remove the adjacent muon's track and energy before 
estimating the amount of hadronic activity nearby. As 
a result, we did not require standard muon isolation in 
this study and instead reduced the overall efficiency by 
a factor of 50% to approximate the loss of signal events 
from modified isolation. 



iD/Cjre V 


le V 


LtlU 


DY+i 


0.15 


0.24 




0.03 


0.08 


tt 


o no 
0.02 


O 1 /I 

U.14 


hi 


< 0.001 


~ 0.03 


T+j 


0.001 


0.002 


IIfJ,+TT 


< 0.001 


< 0.001 


J/^+j 


< 0.001 


< 0.001 


Total 


0.20 


0.49 



TABLE II: Continuum backgrounds for low invariant mass 
muons pairs with missing energy {da / dM^^) for the hP — > 
tt) search at the Tevatron and LHC in units 
of fb/GeV. The backgrounds are given for p^^, _^t, and Ai? 
cuts optimized for a 100 GeV Higgs. 



B. Backgrounds 

There are several backgrounds to this search: Drell- 
Yan muons recoiling against jets, electroweak processes, 
and leptons from hadronic resonances. The Drell-Yan 
background is the most important. The missing energy 
that results from the tau decays is a critical feature in 
discriminating the signal from the background. In ad- 
dition, the fact that the missing energy is in the oppo- 
site direction as the muons reduces the background from 
hadronic semileptonic decays. 

The primary background arises from Drell-Yan 
muons recoiling against a jet. The missing energy is 
either due to mismeasurement of the jet's energy or to 
neutrinos from heavy flavor semi-leptonic decays in the 
jet. In the former instance, the analysis is sensitive to 
how PGS fluctuates jet energies. While PGS does not pa- 
rameterize the jet energy mismeasurement tail correctly, 
the background only needs an 0(30%) fluctuation in the 
energy, which is within the Gaussian response of the de- 
tector. The Drell-Yan background was generated using 
MadGraph/MadEvent , v . 4 . 4 . 16"* I28j and was matched 
up to 3j using an MLM matching scheme. It was then 
showered and hadronized with PYTHIA. Again, the stan- 
dard muon isolation criteria could not be applied and we 
used the same 50% efficiency factor that was used for the 
signal. 

All events are required to have a pair of oppositely- 
signed muons within \ri\ < 2. Each muon must have a 
Pt of at least 10 GeV. A jet veto is placed on all jets, 
except the two hardest. The veto is 15 and 50 GeV 
for the Tevatron and LHC, respectively. Lastly, it is 



^ PYTHIA docs not keep spin correlations in decays. This approx- 
imation does not affect the signal considered here because the 
taus are highly boosted in the direction of aP and any kinematic 
dependence on spin is negligible. As verification of this, TAUOLA 
[26| was used to generate the full spin correlated decays. 



* This version of MadEvent does not apply the xqcut to leptons. 
We thank J. Alwall for altering matrix element-parton shower 
matching for this study. 
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required that the hardest muon and missing energy are 
separated by A(/) > 140° . Table I shows the relative and 
cumulative cut efhciencies for the signal. 

There are three higher-level cuts that further dis- 
tinguish the signal from the background. These cuts are 
optimized as a function of the Higgs mass to maximize 
the significance of the signal. The first is a cut on the 
sum pt of the muons (p^^), and is approximately 



(22) 



The second is a missing energy cut. There is a moderate 
amount of missing energy in the signal events coming 
from the tau decays and this proves to be a very impor- 
tant discriminant from the Standard Model background. 
The missing energy cut is 



^T> (0.2 - 0.25) X rriho. 



(23) 



For the LHC, the requirement is always held above 
30 GeV. The last is a AR cut on the muon pair, which 
depends on both the Higgs and pseudoscalar masses 



Ai?(/i,^) 



> 



(24) 



These cuts depend on the kinematics of the decays and 
the geometry of the events is similar at both the Tevatron 
and LHC. 

Figure [5] shows the invariant mass spectrum for the 
two oppositely-signed muons. The inset shows the signal 
(black) and background (gray) before the pt^^f^, ]^t, and 
AR cuts. After these cuts are placed, the Drell-Yan 
background is mostly eliminated. The muon invariant 
mass reconstructs the mass of the pseudoscalar. 

We used PGS to model the muon invariant mass res- 
olution and used an m„o ±80 MeV to exclude continuum 




5 6 7 8 

Muon Invariant Mass (GeV) 

FIG. 5: Muon invariant mass for 5 fb-i at the LHC before 
(inset) and after the pt^'^ , fir, and AR cuts. The signal, a 
100 GeV Higgs decaying to a pair of 7 GeV pseudoscalars, 
is shown in black and the Drell-Yan background is shown in 
gray. 
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mass for signal with m^o = 
^ at the LHC after p^**, fir, and 



FIG. 6: Total invariant 
100, 150, 200 GeV for 5 fb" 
A_R cuts. The fir is projected in the direction of the hardest 
jet. 



backgrounds. The Drell-Yan background {da / dM^j^^j) is 
0(0.15 fb/GeV) at the Tevatron and 0(0.24 fb/GeV) at 
the LHC. 5 

The other important kinematic handle in this anal- 
ysis is the total invariant mass of the event, which re- 
constructs the mass of the s-channel Higgs boson. The 
total invariant mass of the signal is shown in Fig. [6] af- 
ter all cuts have been applied. The width of the peak is 
narrowed if the missing transverse energy is projected in 
the direction of the jet. We expect that this should be 
the direction of the missing energy because there will be 
boosted neutrinos from the hadronic tau decays. 

In addition to Drell-Yan production, there are sev- 
eral electroweak production mechanisms for muon pairs 
and ]^T- The most important one is W'^W^ production. 
When the vector bosons are in a spin-0 configuration and 
decay leptonically, the muons are nearly coUinear and 
antiparallel to the neutrinos. When the W~ decays to 

the lepton momentum and spin are in the same di- 
rection as the gauge boson. The antineutrino, however, 
is antialigned with the W~ , and thus its momentum is 
antiparallel to that of the muon. The situation is simi- 
lar for the decay, except that the directions of the 
muon and neutrino are reversed. The pL^^~uv back- 
ground was generated with MadGraph and was found to 
be 0(0.03 fb/GeV) at the Tevatron and 0(0.08 fb/GeV) 
at the LHC. 

Top quark production is another important elec- 
troweak background. Using Madgraph to gener- 
ate ^'^ vvbb, we estimate that this background is 
0(0.02 fb/GeV) at the Tevatron and 0(0.14 fb/GeV) at 



The background cross sections we quote in this section are for 
p?^'', and A_R cuts optimized for a ICQ GeV Higgs. 
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the LHC. The top background becomes nearly compara- 
ble to the Drell-Yan background for larger Higgs masses 
due to the weaker p^^, and AR cuts. 

Electroweak production of ^iitt has a production 
cross section on the order of several attobarns when 
requiring low invariant-mass, high pT muons. Conse- 
quently, it is subdominant to the W~^W~ and tt back- 
grounds, with 0(8 X 10"^ fb/GeV) at the Tevatron and 
0{2 X 10-"^ fb/GeV) at the LHC. 

There are several other important backgrounds that 
arise from low-lying hadronic spectroscopy that cannot 
be computed reliably with existing Monte Carlo gener- 
ators. These backgrounds come from (i) semi-leptonic 
decays {b cfj.v), (ii) heavy-flavor quarkonia, and (iii) 
leptonic decays of light mesons. 

Double semileptonic decays {e.g. b c s/d) typ- 
ically give rise to soft leptons in jets but can occasionally 
fluctuate to give hard isolated leptons. It is challeng- 
ing to estimate this background contribution because the 
events are rare and we are statistics- limited. However, 
we have attempted to estimate the relative magnitude 
using PYTHIA. It was found that the total cross section 
for 66 jets to produce two muons is 0(80 fjh) at the LHC. 
Using a power law extrapolation from low pi^'^, it was 
estimated that a pt^^ cut of 40 GeV reduces the cross 
section to 0(10 pb). A missing energy cut of 30 GeV is 
0.6% eflicient and requiring A0(^t', m) > 140° reduces 
the cross section by an additional order of magnitude to 
0(5 fb). Placing a AR cut on the muons and assuming 
that the muon isolation is 10% eflicient, the cross section 
becomes 0(0.3 fb). The invariant mass of the muons in 
these events is distributed over 0(10 GeV), so the fi- 
nal background is approximately 0(0.03 fb/GeV). This 
is likely an overestimate because we have assumed that 
the cuts are uncorrelated. We found that at the Teva- 
tron, semileptonic decays do not produce enough high px 
muon pairs antialigned with the to be an important 
background. 

Ts can decay into muon pairs, but their invariant 
mass is above the range we are interested in. Ts can also 
decay into taus that can subsequently decay into muons 
with a branching fraction of 3%. The invariant mass for 
these muon pairs will be in the region of interest. There 
are, however, two factors that mitigate this background. 
The first is that the muons will be soft unless the T has 
very high p^ ^ 0(60 GeV). The pT spectrum of Ts falls 
off rapidly. At the Tevatron, the differential cross section 
is 0(250 fb) at 20 GeV [29]. A nai've extrapolation 
to 60 GeV would place this background at 0(2 fb). Ad- 
ditionally, the missing energy in these events points in 
the direction of the muon pairs rather than towards the 
recoiling jet and a cut on the angle between the miss- 
ing energy and the muon direction should reduce this 
background by another order of magnitude. Accounting 
for this reduction, as well as the 3% branching fraction 
of the taus to muons, we find that the cross section is 



0(6 X 10-"* fb/GeV) at the Tevatron. We do not expect 
this background to dominate Drell-Yan at the LHC, ei- 
ther; using NNLO predictions for the pt distribution of 
Ts at the LHC |3^, we estimate that this background 
will be 0(2 X 10-3 fb/GeV). 

At the charm threshold, the J/ip and tp{2S) be- 
come important because the tails of distributions aris- 
ing from mismeasurement could spill over into higher 
invariant mass bins. Again, the cross section for the 
decays of these particles drops sharply as a function 
of p^'' and with p':^^^ > 40 GeV, the cross section is 
0(100 fb) [HI]. Because this peak is below the invariant 
mass of interest, only the tail of the distribution 
is a background. The dominant contribution comes ei- 
ther from the Lorentzian tail of the decay width or from 
the non-Gaussian mismeasurement tail. The Lorentzian 
tail suppresses the J/ifj contamination by 0(10^^) for 
m^o between 3.6 and 9 GeV. The Gaussian tail of J/ip 
mismeasurement goes out at least 5a, meaning that the 
contamination should be down by 0(10"^). This gives a 
background cross section smaller than 0(10"^ fb/GeV) 
at the Tevatron and at the LHC. The contributions of 
the '0(25') are subdominant to that of the J/ip [32\. 

Resonances beneath the J/ tp are not a problem be- 
cause they are far enough away from the invariant mass 
window we are interested in. Peaks from fake muons 
may arise from B —f Kt: or similar decays where the 
kaons and pions punch through to the muon chamber. 
These events are typically accompanied by significant 
hadronic activity and tight muon isolation requirements 
(after removing the adjacent muon) will reduce these 
backgrounds of fake muons [33] . Secondly, the I^t from 
in-flight decays is in the direction of the muons, but in 
the signal, it is back-to-back with the muons. Placing a 
cut on the relative angle between the muons and the I^t 
is effective at eliminating these difficult backgrounds. 



C. Expected Sensitivity 

Figure [7] shows the expected 95% exclusion plot at 
the Tevatron and LHC. The contours indicate the cross 
sections for values of (S*)/ sin 2/3; this ratio affects the 
partial width of the Higgs into the pseudoscalars (Eq. 
14). The total projected luminosity for the combined 
data sets at CDF and D0 is 20 fb^^; currently, each 
experiment has ~ 5 fb~^. With 10 fb^^ luminosity, the 
Tevatron will start probing the interesting regime where 
(S) I sin 2/3 = 250 GeV. Once the benchmark luminos- 
ity is reached, the Tevatron will have sensitivity up to 
(S*)/ sin 2/3 = 500 GeV. 

With early data, the LHC has sensitivity to regions 
corresponding to (S")/ sin 2/3 < 250 GeV. The sensitiv- 
ity is weaker for Higgs masses below 100 GeV because 
the backgrounds worsen due to a smaller p^ cut. How- 
ever, combined analyses by CDF and D0 should be able 
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FIG. 7: Expected sensitivity to the Higgs production cross section at the Tevatron (left) and LHC (right) for m^o — 7 GeV. 
The contour hnes indicate the cross sections for several values of (5)/ sin 2/3 (in GeV), which alters Br{h^ a°a°). The 
Standard Model Higgs decay width and (NNLO) gluon fusion production cross sections were obtained from [4]. An e^r of 
0.8% was used for the branching ratio of agao 2^2t. The region beneath the dashed line has been excluded by LEP. 



to probe this region down to a 0(1 pb). By the time 
the LHC reaches a luminosity of 5 fb"^, it will be sensi- 
tive to the most relevant region of parameter space, with 
(S')/sin2/3 < 1 TeV. 

The sensitivity curves depend on the product of the 
pseudoscalar branching ratios into muons and taus, e^,-- 
For Fig. [7] we assumed that the pseudoscalar was 7 GeV, 
which corresponds to e^r — 0.8%. For a Hghter pseu- 
doscalar (e.g., 4 GeV), e^r is nearly double this value. 
In this case, the signal limits can increase by as much as 
a factor of two. To first order in m'^a/'mf^o, the branching 
fraction of the Higgs into the pseudoscalar is independent 
of rriaO and the contour lines in Fig. [7] are unaffected. 
Therefore, if the pseudoscalar is near the tau threshold, 
the experiments are even more sensitive to the Higgs 
production cross section than indicated in the figure. 

IV. CONCLUSION 

We have shown that if the Higgs decays into a 
pair of light pseudoscalars that subsequently decay into 
taus, then the discovery of the Higgs boson is promis- 
ing through the subdominant channel where one pseu- 
doscalar decays to a pair of muons. The Tevatron has 
a chance of discovering this class of models if CDF and 
D(2) perform combined analyses with the full data sets. 
The Tevatron can begin to recover the parameter space 
that LEP missed with their /i" — > 4r search, which was 
prematurely stopped at 86 GeV. Assuming that the only 
new decay mode of the Higgs boson is into a pair of pseu- 
doscalars, the Tevatron is sensitive to m/jo ~ 102 GeV 



with 10 fb"\ and up to m^o ~ 110 GeV with 20 fb"\ 
When the Tevatron covers this ground, their results, 
combined with the direct limits from LEP, will effectively 
establish a lower limit on the Higgs mass regardless of 
the admixture of Higgs decays into light pseudoscalars 
or Standard Model fermions. With a 20 fb~^ cross sec- 
tion, the Tevatron will be sensitive to Higgs bosons up 
to rriho < 150 GeV. 

At the LHC, this search becomes a method of dis- 
covering the Higgs with early data - potentially with 
sub-fb~^ data sets. With an integrated luminosity of 
0(1 fb~^), the LHC will be able to recover the missing 
LEP limits. Eventually, the LHC will be able to push 
this branching ratio down substantially, to the 3% level. 
A discovery or even a limit on such a decay mode will 
be an important step in verifying the field content and 
symmetry structure of the Higgs potential. 
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